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ABSTRACT 

Photolysis of vibrationally excited oxygen produced 
by ultraviolet photolysis of ozone in the upper stratosphere is 
incorporated into the Lawrence Livermore National Laboratory 
2-D zonally averaged chemical -radiative-transport model of the 
troposphere and stratosphere. The importance of this potential 
contributor of odd oxygen to the concentration of ozone is 
evaluated based upon recent information on vibrational distribu- 
tions of excited oxygen and upon preliminary studies of energy 
transfer from the excited oxygen. When the energy transfer 
rate constants of previous work arc assumed, increases in model 
ozone concentrations of up to 40 percent in the upper strato- 
sphere arc found, and the ozone concentrations of the model 
agree with measurements, including data from the Upper Atmo- 
sphere Research Satellite. However, the increase is about 0.4 
percent when the larger energy transfer rate constants suggested 
by more recent experimental work arc applied in the model. 
This indicates the importance of obtaining detailed information 
on vibrationally excited oxygen properties to evaluation of this 
process for stratospheric modeling. 

INTRODUCTION 

Simulations of stratospheric chemistry do not produce 
as large an ozone concentration as is experimentally found in 
the upper stratosphere. This difference has been known for at 
least 15 years (Butler, 1978; WMO, 1985). One of tire possible 
explanations of this lower ozone concentration is that a source 
of odd oxygen is missing from the existing models. Slanger el 
al. (1988) have proposed a possible additional source of odd 
oxygen based on experiment. Photodissocialion of ozone in the 
ultraviolet (200 to 310 nm) produces oxygen molecules in the 
ground electronic state. This oxygen carries a large vibrational 
excitation (Slanger el al ., 1988; Kinugawa el al., 1990). The 
vibrationally excited oxygen (0 2 *) then absorbs a second photon 
by its Schumann-Rungc transition and dissociates to two O 
atoms. This second dissociation can occur with light consider- 
ably to the red of that which phololyzes thermal oxygen. In 
order to incorporate this mechanism into atmospheric models, 
the distribution of 0 2 * with respect to vibrational quantum 
number v and the rale constants for removal of vibrational 
quanta from 0 2 % by collisional quenching must be known. 


Utilizing 0/ distributions from Kinugawa et al. (1990) and 
quenching rate constants for 0 2 * based on Rapp (1965), Toumi 
et al. (1991) found an enhancement of ozone concentration 
ranging from 13 to 60 percent for altitudes of 36 to 58 km in 
their atmospheric model. The 248 nm 0 3 photolysis data of 
Slanger (private communication) indicate a different 0 2 *(v) 
distribution, and preliminary collisional quenching data in the 
same study indicate a considerably higher rate constant than 
was used previously. Photolysis of 0 2 * has been incorporated 
into the LLNL global two-dimensional chemical-radiative- 
transport model, and this study compares the ozone concentra- 
tions calculated using both the earlier collisional quenching and 
the experimental quenching with observed ozone concentrations, 
including recent measurements from the Upper Atmosphere 
Research Satellite Microwave Limb Sounder (UARS MLS). 

EXCITED OXYGEN CHEMISTRY 

The LLNL 2-D model has been discussed previously 
(Johnston et al., 1989). In addition to the standard model, the 
rate of photolysis of 0 3 into the ground state oxygen channel 
(0 2 ( 3 £ f )) is used to calculate 0 2 * processes. As is detailed 
below, the 0 2 4 calculation internally performs three steps for 
each altitude-latitude combination. The nascent distribution of 
0 2 *(v) is first found from 0 3 photolysis rates. Rate constants 
for 0 2 * photolysis arc calculated from absorption cross sections. 
Finally, with the photolysis and energy transfer rate constants, 
the steady stale concentration of 0 2 *(v) is determined. 

When 0 3 is phololyzcd, excess energy is distributed 
among translation (between the centers of mass of 0 2 and O 
atom) and rotation and vibration of the 0 2 fragment. Data 
obtained by Slanger indicate that a bimodal distribution of 
vibrationally excited 0 2 arises when 0 3 is photolyzed at 248 
nm. The time-of-flight study of Kinugawa et al. (1990) shows 
two peaks of O atom with respect to translational energy when 
0 3 is photolyzed at 226 nm. The peak at higher 0 2 vibrational 
quantum number is assigned as partially but not completely due 
to thermal decomposition of 0 3 . We use a double Gaussian fit 
in the relative translational energy of photolysis fragments to 
approximate both the Slanger data and the Kinugawa et al. 
data. Since the vibrational energy in 0 2 ( 5 £ R ) is equal to llie 
photolysis energy minus E,. except for the rotational energy in 
Oj, the vibrational distribution can be calculated. 
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The 0 2 * cross sections arc calculated from wavefunc- 
lions for 0 2 and the Schumann-Rungc transition dipole moment 
function using equations based on those in Levine (1975) and 
Saxon and Slangcr (1991). The wavcfunclions for both the 
ground and excited electronic slates arc calculated from 
potential surfaces using the method of Kulandcr (1988), The 
Rydbcrg-Klcin-Rccs potential surfaces of ground 0 2 ( 3 £ g ) and 
excited 0 2 ( 3 Z„ ) arc from Krupcnic (1972) extrapolated to r = 

1.0 A and to r = 11.0 A. The transition dipole moment 
function D(r) is from Allison et al. (1986). The transition 
strengths produced agree to within 50% with the experimental 
transition strengths of Lewis et al. (1986) and calculated 
transition strengths of Allison et al. (1971). Cross sections 
agree reasonably with the v = 12 cross sections found by Saxon 
and Slangcr (1991). 

The removal of vibrational excitation from 0 2 * 
involves multiple processes which ultimately must be distin- 
guished from one another. The 0 2 * may lose vibrational 
excitation by either transfer to translational and rotational 
energy or transfer to vibrational energy. Either transfer 
mechanism preferentially removes one quantum of vibrational 
energy; however, according to Slangcr 's preliminary observa- 
tions on 0 2 * collisions with nitrogen for v = 18 to 20, processes 
which remove two quanta of vibrational energy from 0 2 * arc 
also significant. Further work at v = 22 finds lhat nitrogen and 
oxygen have approximately the same rate constants for 0 2 * 
quenching, so that v-t transfer must dominate there. A model 
proposed by Rapp and Englander-Golden (1964) and by Rapp 
(1965), as further expanded to multiple-quantum removal by 
Yardlcy (1980), provides the v-v energy transfer rale constants 
used in the 0 2 * model. The v-t contribution chosen for 
compatibility with the Slangcr et al. observations at v = 22 is 

3.0 x 10 15 molcc cm 3 s 1 for all temperatures and for both M 
= 0 2 and M = N 2 ; this is referred to as lire ’Tull k c case” in the 
results. The sensitivity of the 0 2 f processes to the choice of 
energy transfer rate constants must be determined due to the 
preliminary nature of the Slangcr group results. A second 



0 1 0 20 30 40 

0, # Vibrational Quantum Number 


Figure 1. Steady state concentrations of vibrationally excited 
oxygen versus vibrational quantum number for 30° N, 44 km 
on the Autumn Equinox for both zero vibration-to-lranslation 
contribution (o) and the full vibrational energy transfer model 
(□) arc compared with the result of Toumi et al. (1991) for v 
= 15 (a). 


series of 2-D model runs is thus conducted with k^ set to zero 
and M = 0 2 alone; this lower limit to the uncertainty range of 
0 2 * energy transfer rale constants is referred to as the "k^ = 0" 
case below. The rate constants of this case are similar to those 
used by Toumi et al. (1991). 

0 2 * may be either deactivated by collisional energy 
transfer or photolyzed in the Schumann-Rungc transition region. 
The lifetime of 0 2 * is less than one hour, so that the instanta- 
neous equilibrium or steady-stale approximation may be applied 
to each of the species 0 2 *(v). Those equations are solved using 
the fact lhat 0 2 * may not be produced with v greater than its 
dissociation limit of 38. 

MODEL-DERIVED 0 2 * AND 0 3 

Figure 1 illustrates steady state 0 2 * concentrations 
(Equation 4) at noon on the Autumn Equinox (September 22) 
in the LLNL 2-D model with respect to vibrational quantum 
number. In the = 0 model (diamonds), the concentration of 
0 2 * falls off rapidly for v > 25 from values of at least 5 x 10 6 
molcc cm' 5 which remain nearly constant with v. Use of the 
full lq model (squares) reduces the maximum 0 2 * concentration 
to 2 x 10 3 molcc cm' 3 for v = 3, and the concentration decreas- 
es more steeply throughout the range of v, reaching 2 x 
K^molcc cm' 3 at v = 25. The additional quenching of the full 
Iq model severely reduces the calculated 0 2 * concentrations and 
will be shown to reduce the 0 2 * effect upon odd oxygen 
considerably compared with that of Toumi et al. (1991). 
Toumi et al. (1991) provides [0 2 *] only for v = 15 (triangle), 
which is slightly less than is the [0 2 *(15)] for our lq., = 0 case 
in Figure 1. 

The rate of O atom production from 0 2 * also depends 
on the rate of 0 3 photolysis, the solar flux, and the absorption 
cross section of 0 2 *. For both lq models, the rate of 0 2 f 
photolysis in molcc cm' 3 s' 1 (one half of the O atom production 
rate) at 30°N, Autumn Equinox noon is plotted against v in 
Figure 2. The 0 2 * photolysis rate peaks at 2 x 10 5 molcc cm 3 
s ' for v = 20 in the lq. t = 0 model, with considerable contribu- 
tions for v down to 12 and up to 27. However, in the full 
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Figure 2. Rate of photolysis for vibrationally excited oxygen 
versus vibrational quantum number for 30° N, 44 km on the 
Autumn Equinox. Zero vibration-to-translation case: O; full 
vibrational energy transfer model: □; Toumi et al. (1991) for 
v = 15: a. 
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k model, the maximum contribution becomes 1 x 10 3 molcc 
cm J s*’ at v « 12, and the photolysis rate falls off more quickly 
with v to cither direction. The strongest contribution to O atom 
production from 0 2 * comes from v in the range of 10 to 25, and 
the higher v arc most severely affected by use of the full k 
quenching both in terms of concentration and in terms of 
photolysis rale. 

Figure 3 illustrates the ratio of the 0 2 * dissociation 
rate to the thermal 0 2 dissociation rate at the Autumn Equinox. 
The altitude shown is the log-pressure scaled altitude used in 
the 2-D model. Altitudes less than 30 km are omitted for this 
figure and Figure 4 since neither process is significant. In the 
k,., = 0 case (Figure 3 A), the 0 2 * process contributes 10 percent 
as much odd oxygen as docs thermal 0 2 photodissociation 
around 35 km for all latitudes except the Antarctic region. The 
odd oxygen addition from 0 2 * increases with increasing altitude 
up to 55 km. Peak values of the ratio arc greater than 50 
percent at latitudes of 45 °S and 60 °N and 55 km altitude, 
Toumi et ai (1991) Table 2 shows an odd oxygen production 
ratio at 30 °N on the Autumn Equinox which increases mono- 
tonically from 13 percent at 36 km to 60 percent at 58 km, but 
with a leveling off of the rate of increase with altitude at the 
higher altitudes. Some quantitative differences in the produc- 
tion ratio exist between our model and that of Toumi et ai, but 
the qualitative features arc mostly the same at 30 °N. 



Figure 3. Ratio (percent) of rate of photolysis (molcc cm 3 s' 1 ) 
of 0 2 * to that of thermal 0 2 for (A) k t = 0 case; (B) full k c 
case. 



Figure 4. Comparison of LLNL 2-D model results with 
measurements of ozone concentration on the Autumnal Equinox 
at 30° N. * - UARS MLS (lines represent lo uncertainty); + 

- MAP recommended values (10 percent typical uncertainty); a 

- LLNL 2-D ambient atmosphere (no 0 2 *); O - LLNL 2-D 
model, k t = 0.0 cm 3 molcc 1 s \ 0 2 v-v transfer only; □ - 
LLNL 2-D model, k t = 3.0 x 10 15 cm 3 molec 1 s ', N 2 and 0 2 
colliders. 

The full k case (Figure 3B) reduces die odd oxygen 
percentage ratio from 0 2 * compared with dicrmal 0 2 by two 
orders of magnitude for all latitudes and altitudes. Increased 
collisional quenching forces the maximum 0 2 * production 
upward (lower pressure), so that the maximum of 0.4 percent 
is located at 58 km, 60 °S or 60 °N; the upper limit of die 
model is at 58 km, so that this is not necessarily die overall 
maximum, but contributions for higher altitudes should not be 
markedly larger due to decreasing concentrations of () 3 . The 
contribution of 0 2 * below 58 km is sharply reduced at all 
latitudes toward a minimum of less than 0.1 percent about 40 
km. 0 3 concentrations for an ambient atmosphere in our 2-D 
model also typically change by 0.1 percent from year to year. 

Addition of approximately half again as much O atom 
from 0/ as is produced by dicrmal 0 2 photolysis in the k- t = 
0 case to the upper stratosphere should produce a large increase 
in 0 3 concentration. Figure 4 compares model 0 3 concentra- 
tions on the Autumn Equinox at 30°N as a function of altitude 
with preliminary results from UARS MLS and the recommend- 
ed zonal mean values of the Middle Atmosphere Program 
(Kcaling el ai., 1989). When k t = 0 is used (diamonds), the 
0 3 concentration is markedly greater for all altitudes above 32 
km than for die ambient atmosphere. The resulting concen- 
trations arc widiin the uncertainly limits of the preliminary 
measurements of die UARS MLS and arc comparable to die 
MAP recommended values for September at 30 °N (crosses) for 
log-pressure altitudes down to 44 km. The MLS measurements 
arc an average over 28 °N to 32 °N for September 21, 1991 for 
the 205 GHz 0 3 band. Absolute accuracy of MLS is still under 
investigation, but the 0 3 concentrations observed at these 
altitudes arc not expected to vary significandy. The change 
between the LLNL 2-D model without 0 2 * processes (triangles) 
and for the 0 2 * full k casc (squares) is insignificant on the 
scale of Figure 4. 
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CONCLUSIONS 

Odd oxygen production from the photolysis of 
vibrationally excited oxygen (0 2 *) produced by ultraviolet 
photolysis of ozone (0 3 ) in the upper stratosphere has been 
incorporated in the LLNL2-D global stratosphere- trophospherc 
model. When quenching of O a * only by v-v collisional energy 
transfer to 0 2 is considered, odd oxygen production rales 
throughout the upper stratosphere (>30 1cm) increase by up to 
50 percent, with a corresponding increase in 0 3 concentration, 
in qualitative agreement with the previous work by Toumi et al. 
(1991). Inclusion of an assumed v-t collisional energy transfer 
rate constant and of N 2 as a collider gas, indicated by recent 
preliminary experimental work, reduces the odd oxygen 
contribution from 0 2 * processes to insignificance compared with 
thermal 0 3 photolysis. The sensitivity of 0 2 * processes in the 
stratosphere to detailed energy transfer rate constants is well 
demonstrated here; unfortunately, these constants arc not yet 
well determined in the laboratory. Thus, this study shows the 
importance of obtaining detailed experimental information on 
the input parameters of O/, particularly the rate constant for 
collisional energy transfer as a function of v, in order to judge 
the potential importance of 0 2 t processes to stratospheric 0 3 
concentrations. 

ACKNOWLEDGMENTS 

This work was performed under the auspices of the U. 
S. Department of Energy by die Lawrence Livermore National 
Laboratory under contract No. W-7405-Eng-48 and was 
supported in part by the U. S. Department of Energy (DOE) 
Office of Health and Environmental Research’s Environmental 
Sciences Division and by die National Aeronautics and Space 
Administration Upper Atmosphere Research Satellite program. 
KP also appreciates appointment to the Global Change Distin- 
guished Postdoctoral Fellowships sponsored in part by the U. S. 
DOE, Office of Health and Environmental Research, and 
administered by Oak Ridge Associated Universities. We thank 
R. P. Saxon and D. L. Hucstis at SRI International, Inc. for 
many helpful discussions concerning this study, and K. Kulan- 
dcr and B. P. Lcngsficld of LLNL for use of and guidance in 
their vibrational wavcfunction code. 

REFERENCES 

Allison, A. C., A. Dalgamo, and N. W. Pasachoff, 
1971: Absorption by Vibrationally Excited Molecular Oxygen 
in the Schumann- Run gc Coplinuum Planet. Space Sci. 19, 
1463-73. 

Allison, A. C., S. L. Guberman, and A. Dalgamo, 
1986: A Model of the Schumann-Rungc Continuum of 0 2 J . 
Geophys. Res. 91, 10,193-8. 

Butler, D. M., 1978: The Uncertainty in Ozone 

Calculations by a Stratospheric Photochemistry Model Geo - 
phys. Res. Lett. 5, 769-72. 

, Johnston, H. S., D. E. Kinnison, and D. J. Wucbblcs, 

1989: Nitrogen Oxides from High-Altitude Aircraft: An 

Update of Potential Effects on Ozone J. Geophys. Res. 94, 
16,351-63. 

Keating. G. M., M. C. Pitts, and D. F. Young, 1989: 
Ozone Reference Models for lire Middle Atmosphere (New 


CIRA) in Middle Atmosphere Program: Handbook for MAP 
31, G. M. Keating cd., pp. 1-36. 

Kinugawa, T., Sato, T., Arikawa, T., Matsumi, Y., 
and Kawasaki, M., 1990: Formation of 0( 3 Pj) Photo fragments 
from the Hartley Band Photodissociation of Ozone at 226 nm 
J . Chem. Phys. 93, 3289-94. 

Krupcnic, P. H., 1972: The Spectrum of Molecular 
Oxygen J. Phys. Chem. Ref. Data 1, 423-534. 

Kulandcr, K. C., 1988: Spectroscopy of SO Chem. 
Phys. Lett. 149. 392-6. 

Levine, I. N., 1975: Molecular Spectroscopy (John 
Wiley and Sons, New York), pp. 306-8. 

Lewis, B. R., L. Berzins, and J. H. Carver, 1986: 
Oscillator Strengths for the Schumann-Runge Bands of ,& 0 2 7. 
Quant. Spectrosc. Radiat. Transfer 36, 209-32. 

Rapp, D„ and P. Englander-Golden (1964): Resonant 
and Near-Resonant Vibrational- Vibrational Energy Transfer 
between Molecules in Collisions J. Chem. Phys. 40, 573-5. 

Rapp, D., 1965: Interchange of Vibrational Energy 
between Molecules in Collisions J. Chem. Phys. 43, 316-7. 

Saxon, R. P., and T. G. Slanger, 1991: Relative 

Contributions of Discrete and Continuum Absorption to the 
Photodissociation of Vibrationally Excited 0 2 (v = 12-20) J. 
Geophys. Res. 96, 17291-5. 

Slanger, T. G., L. E. Jusinski, G. Black, and G. E. 
Gadd, 1988: A New Laboratory Source of Ozone and Its 

Potential Atmospheric Implications Science 241, 945-50. 

Toumi, R., B. J. Kerridge, and J. A. Pyle, 1991: 
Highly Vibrationally Excited Oxygen as a Potential Source of 
Ozone in the UppeT Stratosphere and Mesosphere Nature 351. 
217-9. 

WMO, 1985: Comparison of Calculated and Ob- 

served Ozone Profiles in Atmospheric Ozone 1985: Assessment 
of our Under st arvling of the Processes Controlling its Present 
Distribution and Change (World Meteorological Organization 
Report Number 16), vol. 2, section 8.2, pp. 420-9. 

Yardley, J. T„ 1980: Introduction to Molecular 

Energy Transfer (Academic Press, New York). 

Note added in proof After this work was presented, R. Toumi 
has published "An Evaluation of Aulocatalytic Ozone Produc- 
tion from Vibrationally Excited Oxygen in the Middle Atmo- 
sphere" in/. Atmos. Chem, 15, 69-77. That report incorporated 
past experimental quenching results and found ozone enhance- 
ments ranging from 0.4 to 1 1 percent for 40-60 km altitude, 
depending on various models of the nascent 0 2 * distribution 
from 0 3 photolysis. 
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